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The effects of the spin-orbit interaction on the tunneling magnetoresistance of magnetic tun- 
nel junctions are investigated. A model in which the experimentally observed two-fold symmetry 
of the anisotropic tunneling magnetoresistance (TAMR) originates from the interference between 
Dresselhaus and Bychkov-Rashba spin-orbit couplings is formulated. Bias induced changes of the 
Bychkov-Rashba spin-orbit coupling strength can result in an inversion of the TAMR. The the- 
oretical calculations are in good agreement with the TAMR experimentally observed in epitaxial 
Fe/GaAs/Au tunnel junctions. 
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The tunneling magnetoresistance (TMR) effect refers 
to ferromagnet/insulator/ferromagnet heterojunctions, 
in which the magnetoresistance exhibits a strong 
dependence of the relative magnetization directions 
in the different ferromagnetic layers and their spin 
polarizations ,ii2i^ Because of this peculiarly strong asym- 
metric behavior of the magnetoresistance, TMR devices 
find multiple uses ranging from magnetic sensors to mag- 
netic random access memory applications^^ 

It came as a surprise that the tunneling magnetore- 
sistance may strongly depend on the absolute orienta- 
tion of the in-plane magnetization directions with respect 
to a fixed crystallographic axis, as experimentally ob- 
served in Refs. The phenomenon was termed tun- 
neling anisotropic magnetoresistance (TAMR)M Even 
more intriguing is the observation of the TAMR effect 
in heterojunctions such as (Ga,Mn)As/A10x/Ar£ and 
Fe/GaAs/Au^ sandwiches, where only one of the lay- 
ers is magnetic, and for which the TMR effect is ab- 
sent. It has been recognized that the origin of the 
TAMR effect is related to the spin-orbit interaction 
(SOI) However, the nature and details of the 
underlying mechanism producing the TAMR remains a 
puzzle. In fact, it has become clear that the responsible 
mechanisms for the TAMR can be different in different 
systems. The (Ga,Mn) As/alumina/ Au heterojunctions^ 
and in (Ga,Mn)As nanoconstrictionsii has been associ- 
ated with the anisotropic density of states in the ferro- 
magnet (Ga.Mn)As and was theoretically modelled by 
introducing strain effects. First-principle calculations for 
the case of an Fe(OOl) surface have recently demonstrated 
the appearance of the TAMR effect due to the shifting 
of the resonant surface band via Rashba SOI when the 
magnetization direction changes^ Furthermore, it has 
recently been observed that the symmetry axis of the 
two- fold symmetry of the TAMR in Fe/GaAs/Au hetero- 
junctions can be flipped by changing the bias voltage^ 

Here we formulate the model proposed to explain 
the experimental results of Ref. |9j, in which the two- 
fold symmetry of the TAMR observed in epitaxial fer- 
romagnct/scmiconductor/normal metal junctions origi- 
nates from the interface-induced Civ symmetry of the 



SOI arising from the interference of Dresselhaus and 
Bychkov-Rashba spin-orbit couplings. Such interference 
effects has already been investigated in lateral transport 
in 2D electrons systems^ The symmetry, which is im- 
printed in the tunneling probability becomes apparent in 
the contact with a magnetic moment. 

Consider a ferromagnet/semiconductor/normal-metal 
tunnel heterojunction. The semiconductor is assumed to 
lack bulk inversion symmetry (zinc-blende semiconduc- 
tors are typical examples). The bulk inversion asymme- 
try of the semiconductor together with the structure in- 
version asymmetry of the heterojunction give rise to the 
Dresselhau o 14 i 15 i 16 and Bychkov-Rashba^ SOIs, respec- 
tively. The interference of these two spin-orbit couplings 
leads to a net, anisotropic SOI with a Ci v symmetry 
which is transferred to the tunneling magnetoresistance 
when the electrons pass through the semiconductor bar- 
rier. The model Hamiltonian describing the tunneling 
across the heterojunction reads 



Here 



H = Hq + Hz + Hbr + Hi 
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with V(z) the conduction band profile defining the po- 
tential barrier along the growth direction (z = [001]) of 
the heterostructure. The electron effective mass m(z) 
is assumed to be m = m c in the central (semiconductor) 
region and m — mi = m r mo (here mo is the bare elec- 
tron mass) in the left (ferromagnetic) and right (normal 
metal) layers. 

The spin splitting due to the exchange field in the fer- 
romagnetic layer is given by 
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Here A(z) represents the exchange energy, er is a vector 
whose components are the Pauli matrices, and n is a unit 
vector defining the spin quantization axis determined by 
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the in-plane magnetization direction in the ferromagnct. 
The Zeeman splitting in the semiconductor and normal 
metal can be neglected. 

The Dresselhaus SOI can be written a3 I5 i i6 1 iLi8 1 i9 
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where x and y correspond to the [100] and [010] direc- 
tions, respectively. The Dresselhaus parameter 7(2) has 
a finite value 7 in the semiconductor region, where the 
inversion bulk inversion asymmetry is present, and van- 
ishes elsewhere. Note that because of the step- like spatial 
dependence of j(z), the Dresselhaus SOI [Eq. (|4])] implic- 
itly includes the interface and bulk contributions, 1 - 5 

The Bychkov-Rashba SOI due to the interface inver- 
sion asymmetry is incorporated in the model through the 
termM 



Hbr = 7^] ai(a x p y - a v p x )5(z - Zi), 



(5) 



where, ai (a r ) denotes the SOI strength at the left (right) 
interface zi = (z r = d). The Bychkov-Rashba SOI 
contribution inside the semiconductor can be neglected 
here. 

Assuming that the in-plane wave vector k|| is con- 
served throughout the heterostructure, one can decou- 
ple the motion along the growth direction (z) from the 
other spatial degrees of freedom. The z component of the 
scattering states in the left (ferromagnetic) region [eigen- 
states of the Hamiltonian |T])] with eigenenergy E can be 
written as 
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where z < 0, x<j represents a spinor corresponding to a 
spin parallel (a =|) or antiparallcl (a =1) to the magne- 
tization direction defined by the vector n, and k a is the 
corresponding z component of the wave vector in the left 
region. In the central (semiconductor) region (0 < z < d) 
we hav o 17 ' 19 
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where q± = (1 T Sm^^n/ft 2 ) -1 / 2 ^ (with q being the 
z component of the wave vector in the barrier in the 
absence of SOI) and £± are spinors corresponding to spins 
parallel (+) and antiparallcl (— ) to the direction k|| x z, 
which is the quantization direction in the barrier. In the 
right (normal metal) region [z > 0) the scattering states 
read 



(r) 



t e *«*<>-*> Xa + t^e^-'^x-a, (8) 



where K a is the corresponding z component of the wave 
vector in the right region. The expansion coefficients in 



Eqs. d6j) - (|8]) can be found by applying standard match- 
ing conditions at each interface i 19 ' 21 Once the wave func- 
tion is determined, the particle transmissivity can be cal- 
culated from the relation 

T„(E,k\\) = Re[n a \t«, a \ 2 + K _ a \t a ^ a \ 2 }. (9) 

The current flowing along the hetcrojunction then is 

1 = E J - ( 10 ) 

where 

Ia = jdrn J dEd2k \M E ^\\)lh( E ) - fr(E)}, (11) 

and fi(E) and f r (E) are the electron Fermi-Dirac dis- 
tributions with chemical potentials [i\ and \i r in the left 
and right leads, respectively. 

The TAMR refers to the changes of the tunneling mag- 
netoresistancc (R) when varying the magnetization direc- 
tion n of the magnetic layer with respect to a fixed axis. 
Here we assume the magnctorcsistance R[uo], measured 
when n points in the [110] crystallographic direction as 
a reference. The TAMR is then given by 

TAMR [110]W = fl(g) „-^ 10 U J tuoi-^) ; (12) 
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where 9 is the angle between the magnetization direction 
n = (cos 8, sin 8, 0) and the [100] crystallographic axis. 
We also find it useful to define the tunneling anisotropic 
spin polarization (TASP) as 



TASP [11O] (0) 
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The TASP measures the changes in the tunneling spin 
polarizatio n 2 ' 4 ' 17 P = (if — Ii)/I (which is a measurable 
quantity^ accounting for the polarization efficiency of the 
transmission) when rotating the in-plane magnetization 
in the ferromagnet. 

For a concrete demonstration of the proposed theoret- 
ical model we performed zero-temperature calculations 
of the TAMR in an epitaxial Fe/GaAs/Au heterojunc- 
tion similar to that used in the experimental observations 
reported in Ref. H|. We use the value m c = 0.067 m 
for the electron effective mass in the central (GaAs) re- 
gion. The barrier width and high (measured from the 
Fermi energy) are, respectively, d = 80 A and V c = 
0.75 eV, corresponding to the experimental samples in 
Ref. 0. For the Fe layer a Stoner model with the ma- 
jority and minority spin channels having Fermi momenta 
/cft = 1-05 x 10 8 cm" 1 and k Fl = 0.44 x 10 s cm" 1 ^ 
respectively, is assumed. The Fermi momentum in Au is 
taken asKi? = 1.2xl0 8 cm -1 .— We consider the case of 
relatively weak magnetic fields (specifically, B = 0.5 T). 
At high magnetic fields, say, several Tesla, our model is 
invalid as it does not include cyclotron effects relevant 
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FIG. 1: Values of the ratio -R[iio]/-R[iio] as a function of the 
interface Bychkov-Rashba parameters ai and a r - 



when the cyclotron radius becomes comparable to the 
barrier width. 

The Dresselhaus spin-orbit parameter in GaAs is 7 = 
24 cV k 3 M^JA The values of the Bychkov-Rashba pa- 
rameters ai, a r [see Eq. (JS])] are not know for metal- 
semiconductor interfaces. Due to the complexity of the 
problem, a theoretical estimation of such parameters 
requires first principle calculations including the band 
structure details of the involved materials, which is be- 
yond the scope of the present paper. Here we assume 
ai and a r as phcnomenological parameters. In order to 
investigate how does the degree of anisotropy depend on 
these two parameters we performed calculations of the 
ratio iJfiio] /R[uo] (which is a measure of the degree of 
anisotropy^) as a function of ai and a r . The results are 
shown in Fig. [TJ where one can appreciate that the size 
of this ratio (and, consequently, of the TAMR) is dom- 
inated by ai. Then, since the values of the TAMR arc 
not very sensitive to the changes of a r we can set this pa- 
rameter, without loss of generality, to zero. This leaves 
ai as a single fitting parameter when comparing to ex- 
periment. Such a comparison is shown in Fig. Ufa) for 
different values of the bias voltage. The agreement be- 
tween theory and experiment is very satisfactory. The 
values of the phcnomenological parameter ai are deter- 
mined by fitting the theory to the experimental value of 
the ratio i?[iio]/-R[iio] an( i this is enough for our theoreti- 
cal model to reproduce the complete angular dependence 
of the TAMR, i.e., the proposed model is quite robust. 
Assuming that the interface Bychkov-Rashba parameters 
are voltage dependent^ (unlike 7, which is a material 
parameter) we perform the same fitting procedure for 
the available experimental data corresponding to differ- 
ent bias voltages and extract the bias dependence of ai 
[see the inset in Fig. Eta)], from which we estimate the 

o 2 

value ai f=a 23 meV A at zero bias. This value is compa- 
rable to the corresponding value of the interface Bychkov- 
Rashba parameter ai f=a 27 meV A obtained from a five 
level k- p model^ for an InAs/GaAs interface. Selecting 
InAs/GaAs for comparison with Fe/GaAs we only wish 
to show that our fitted values have reasonable magnitude, 
not differing too much from known values in scmiconduc- 




FIG. 2: (a) Angular dependence of the TAMR in an 
Fe/GaAs/ Au tunnel heterojunction for different values of the 
bias voltage Vuas- Solid lines corresponds to our theoret- 
ical results while symbols represent the experimental data 
(conveniently mirrored) as deduced from Ref. U The values 
of the phenomenological parameter on have been determined 
by fitting the theory to the experimental values of the ratio 
R\lio] / R[iio] for each value of Vuas- The extracted bias de- 
pendence of ai is shown in the inset, (b) Angular dependence 
of the TAMR for different barrier widths. 



tor interfaces. Interestingly ai in our system changes sign 
at a bias slightly below 50 mV. This bias induced change 
of the interface Bychkov-Rashba parameter results in an 
inversion of the TAMR [see Fig. [Ha)] . Similar behavior 
is reported by ab initio calculations on Fe surfaces, where 
only Bychkov-Rashba SOI is present^ 

The robustness of our model can be understood from 
the following simplified picture of the TAMR effect. The 
SOI term Hso = Hq + Hbr can be written [see Eqs. (U]) 
and ([5])] as a Zeeman-like term Hso ~ B e / / • cr with the 
effective magnetic field 

B e //(k||) = (ai5{z)k y - jk x dl, -ai5{z)k x + jk v dl,0), 

(14) 

where, for the sake of qualitative argument we neglect 
the interface Dresselhaus contributions. Performing the 
average of B e // over the unperturbed (in the absence 
of SOI) states of the system one obtains the following 
general form of the averaged effective spin-orbit magnetic 
field 

w(k||) = (aik y - jk x ,-aik x +jk y ,0), (15) 

where ai ~ aif a (k^) and 7 = 7/7(^11), with f a (k\\) and 
/ 7 (fcy) being real functions of k\\ = |k|||. The effective 
field w(k||) becomes anisotropic in the ky-space with a 
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C2v symmetry when both ai and 7 have finite values. 
It characterizes the amount of k|| -dependent precession 
of the electron spin during the tunneling. For a given 
k|| there are only two preferential directions in the sys- 
tem, defined by n and w(k||). Therefore, the anisotropy 
of a scalar quantity such as the total transmissivity 
T(S,k||) = T T (S,k||) + Ti(E,kn) can be obtained as a 
perturbative expansion in powers of n • w(k||), since the 
SOI is much smaller than the other relevant energy scales 
in the system. The total transmissivity is then given, 
up to second order in the anisotropy, by T(£ , ,kii) as 
TW^fci^+TW^fci^n-wCk^+TW^fci^n-wCkn)] 2 . 
Averaging over the in-plane momenta to get the full cur- 
rent, one obtains for small bias voltages, 

/ = (r( )(^F,fc||)) + (r( 2 )(^,fc||)[n-w(k||)] 2 )n ws , 

(16) 

where (...) represents average over k||. Note that the first 
order term vanishes after average over k|| since w(k||) = 
-w(-k||). Taking into account Eqs. |[T2"]). ([T5|). and fT6]) 
one obtains the following approximate expression for the 
TAMR 

(&aTWkf,)[sm(29) - 1] 
TAMR [110] » jL cx a,7[8in(2fl)-l], 

(17) 

where the arguments of the expansion coefficients 
and T^ 2 ) have been omitted for brevity. The angular 
dependence in Eq. (|17p is consistent with that found ex- 
perimentally, as well as that obtained from the full the- 
oretical calculations [see Fig. Ufa)]. One can clearly see 
from Eq. (jTTJ) that bias-induced changes of the sign of 
the Bychkov-Rashba parameter u\ lead to an inversion 
of the TAMR. When a n = 0, the two-fold TAMR is 
suppressed. To put in words, the Bychkov-Rashba (or 
Drcsselhaus) term by itself cannot explain the observed 
C2V symmetry. The suppression of the TAMR is approx- 
imately realized in Fig. Ufa) for the case of a bias voltage 
of 50 mV. 

The above model neglects the contribution of the spin- 
orbit-induced symmetries of the involved bulk structures. 
Say, Fe exhibits a four-fold anisotropy, which should be 
reflected in the tunneling density of states. The fact that 



this is not seen in the experiment suggests that this effect 
is smaller than the two-fold symmetry considered in our 
model. 

Another system parameter that can influence the size 
of the TAMR is the width of the barrier. The angular de- 
pendence of the TAMR for the case of Vuas = — 90 meV 
is displayed in Fig. [D(b) for different values of the barrier 
width d. As clearly seen in Fig. [2Kb), our model predicts 
an increase of the TAMR when increasing the width of 
the barrier. 

Finally, we show the angular dependence of the TASP 
[see Eq. (|13p ] in Fig.[3]for different values of the bias volt- 
age. The anisotropy of the tunneling spin polarization 
indicates that the amount of transmitted and reflected 
spin at the interfaces depends on the magnetization di- 
rection in the Fc layer, resulting in an anisotropic spin 
[110] [-100] [-110] 
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FIG. 3: Angular dependence of the TASP for different values 
of the bias voltage. 

local density of states at the Fermi surface 2 ^ and showing 
spin- valve-lake characteristics. 

In summary, we have formulated a theoretical model in 
which the two-fold symmetry of both the TAMR and the 
TASP in epitaxial ferromagnet /semiconductor/normal- 
metal heterojunctions originates from the interplay be- 
tween the Dresselhaus and Bychkov-Rashba SOIs. Our 
theoretical results for epitaxial Fe/GaAs/Au heterojunc- 
tions are in good agreement with the available experi- 
mental data. 
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